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Abstract 

We study the limit of high activation energy of a special Fokker-Planck equation, known 
as Kramers-Smoluchowski (K-S) equation. This equation governs the time evolution of the 
probability density of a particle performing a Brownian motion under the influence of a 
chemical potential H/e. We choose H having two wells corresponding to two chemical states 
A and B. We prove that after a suitable rescaling the solution to (K-S) converges, in the limit 
of high activation energy (e — > 0), to the solution of a simple system modeling the diffusion 
of A and B, and the reaction A^ B. 

The aim of this paper is to give a rigorous proof of Kramer's formal derivation and to 
embed chemical reactions and diffusion processes in a common variational framework which 
allows to derive the former as a singular limit of the latter, thus establishing a connection 
between two worlds often regarded as separate. 

The singular limit is analysed by means of Gamma-convergence in the space of finite Borel 
measures endowed with the weak-* topology. 

Key words and phrases: unification, scale-bridging, upscaling, high-energy limit, acti- 
vation energy, Dirichlet forms, Mosco-convergence, variational evolution equations 
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1 Introduction 

1.1 Chemical reaction as a diffusion process 

In a seminal paper in 1940, Hendrik Anthony Kramers described a number of approaches to 
the problem of calculating chemical reaction rates [T^ . One of the limit cases in this paper is 
equivalent to the motion of a Brownian particle in a (chemical) potential landscape. In this 
description a reaction event is the escape of the particle from one energy well into another. 

This description is interesting for a number of reasons. It provides a connection between 
two processes, diffusion and reaction, which are often — especially at the macroscopic level — 
viewed as completely separate. It also provides a link between a macroscopic effect — chemical 
reaction — and a more microscopic, underlying motion, and in doing so, it highlights the fact 
that diffusion and reaction ultimately spring from the same underlying motion. It finally also 
allows for explicit calculation of reaction rates in terms of properties of the energy landscape. 

In this paper we contribute to this discussion by studying the limit process of high acti- 
vation energy in the unimolecular reaction A ^ B. As a first contribution, this provides a 
rigorous proof of the result that Kramers had derived formally. At the same time we extend 
his result to a Brownian motion in the product space spanned by both the chemical variable 
of Kramers and the variables corresponding to position in space, resulting in a limit system 
that models not only chemical reaction but also spatial diffusion — a simple reaction-diffusion 
system. 
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With this paper we have two aims. The first is to clarify the mathematical — rigorous — 
aspects of the formal results of |12) . and extend them to include spatial diffusion, and in this 
way to contribute to the upscaling of microscopic systems. The second is to make a first step 
in the construction of a variational framework that can describe the combination of general 
diffusive and chemically reactive processes. From this point of view it would be interesting, 
for example, to place the limit system in the context of Wasserstein gradient flows (see also 



Section 1.10 1. Initiated by the work of Otto [111 I15| and extended into many directions since, 
this framework provides an appealing variational structure for very general diffusion processes, 
but chemical reactions have so far resisted representation in the Wasserstein framework. 

In this paper we only treat the simple equation A ^ B, but we plan to extend the approach 
to other systems in the future (see also [14)'). 

1.2 The setup: enthalpy 

We consider the unimolecular reaction A B. In chemical terms the A and B particles are 
two forms of the same molecule, such that the molecule can change from one form into the 
other. A typical example is a molecule with spatial asymmetry, which might exist in two 
distinct, mirror-image spatial configurations; another example is that of enzymes, for which 
the various spatial configurations also have different biological functions. 

Remark. Classical, continuum-level modelling of the system of A and B particles that diffuse 
and react (see e.g. [SjISj) leads to the set of differential equations, where we write A and B 
for the concentrations of A and B particles: 

dtA~ DAA = k{B - A) (la) 
dtB - DAB = k{A- B). (lb) 

(See Section [l.lO| for the equal reaction rates). This system will arise as the upscaling limit 
(see Theorem yj) of the system that we now develop in detail. 

We next assume that the observed forms A and B correspond to the wells of an appro- 
priate energy function. Since it is common in the chemical literature to denote by 'enthalpy 
difference' the release or uptake of heat as a particle A is converted into a particle B, we shall 
adopt the same language and consider the A and B states to correspond to the wells of an 
enthalpy function H. 

While the domain of definition of H should be high-dimensional, corresponding to the 
many degrees of freedom of the atoms of the molecule, we will here make the standard 
reduction to a one-dimensional dependence. The variable ^ is assumed to parametrize an 
imaginary 'optimal path' connecting the states A to B, such that ^ = — 1 corresponds to 
A and ^ = 1 to _B. Such a path should pass through the 'mountain pass', the point which 
separates the basins of attraction of A and B, and we arbitrarily choose that mountain pass 
to be at 5 = 0, with H{0) = 1. We also restrict ^ to the interval [—1, 1], and we assume for 
simplicity that the wells are at equal depth, which we choose to be zero. A typical example 
of the function H is showed in Figure [l] 
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Figure 1: A typical function H 

Specifically we make the following assumptions about H: H £ C°°([— 1, 1]), H is even in 
5, maximal at ^ = with value 1, and minimal at ^ = ±1 with value 0; -ff(C) > for any 
— 1 < ^ < 1; _ff'(±l^) = 0. The assumption of equal depth for the two wells corresponds to an 
assumption about the rate constants of the two reactions; we comment on this in Section [l.lO| 
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1.3 Diffusion in the chemical landscape 



This newly introduced 'chemical variable' ^ should be interpreted as an internal degree of 
freedom of the particle, associated with internal changes in configuration. In the case of two 
alternative states of a molecule, ^ parametrizes all the intermediate states along a connecting 
path. 

In this view the total state of a particle consists of this chemical state ^ together with 
the spatial position of the particle, represented by a d-dimensional spatial variable a; in a 
Lipschitz, bounded, and open domain C K'', so that the full state space for the particle is 
the closure D of 

_D := S7 X (—1, 1) with variables (x,^). 

Taking a probabilistic point of view, and following Kramers, the motion of the particle will 
be described in terms of its probability density p £ ^{D), in the sense that for Borel sets 
X C n and H C [—1,1] the number p{X x H) is the probability of finding the particle at a 
position X (z X and with a 'chemical state' ^ € H. 

The particle is assumed to perform a Brownian motion in D, under the influence of the 
potential landscape described by H. This assumption corresponds to the 'large- friction limit' 
discussed by Kramers. The time evolution of the probability distribution p then is given by 
the Kramers-Smoluchowski equation 

dtp-A^p~rdi,{di:p + pd^H)=0 in X (0,oo)), (2) 

with Neumann boundary conditions on the lateral boundary dD. The coefficient r > is 
introduced to parametrize the difference in scales for x and ^: since a: is a rescaled physical 
distance, and ^ is a rescaled 'chemical' distance, the units of length in the two variables are 
different, and the parameter r can be interpreted as the factor that converts between the two 
scales. Below we shall make an explicit choice for r. 

1.4 The limit of high activation energy 

In the setup as described above, there is a continuum of states (i.e. (—1, 1)) connecting the A 
state to the B state, and a statement of the type 'the particle is in the A state' is therefore not 
well defined. In order to make a connection with the macroscopic description 'A B', which 
presupposes a clear distinction between the two states, we take the limit of high activation 
energy, as follows. 

We rescale the enthalpy H with a small parameter e, to make it H{^)/e. (This is called 
'high activation energy' since maxj H{^)/e = 1/e is the height of the mountain that a particle 
has to climb in order to change states). 

This rescaling has various effects on the behaviour of solutions p of ([2| . To illustrate one 
effect, let us consider the invariant measure 7^, the unique stationary solution in 3^{D) of ([2]): 

7. = An® 7., ^"-^^^^'in' 7. = ^r'e-^/"£'|j_^ ^ (3) 

(where £^,£'' are the 1- and d-dimensional Lebesgue measures). The constant is fixed by 
the requirement that 7e(-D) = 7£([— 1, 1]) = 1. 
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Figure 2: The density % 

Since H is strictly positive at any —1 < ^ < 1, the exponential exp(— //(^)/e) vanishes at 
all ^ except for ^ = ±1; therefore the measure 7^ concentrates on the lines ^ = — 1 and ^ = 1, 
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and converges weakly-* as e — > to the limit measure 7 given by 

7 = An®7, ^ ■= + Si). 



(4) 



Here weak-* convergence is to be interpreted in the duality with continuous functions in D 
(thus considering ^{D) as a weakly-* closed convex subset of the space ^{D) = (CP{D)) 
of signed Borel measures with finite total variation) i.e. 

Mm f_^{x,Od'y,^ f_(^{x, djix, = I f [H^, ~l)+Hx,l)) dXn{x), for any <^ £ C"(:D). 
'^^^ J D J D ^ Jn 

We should interpret the behaviour of 7^ as follows. In the limit e ^ 0, the deep wells at 
^ = ±1 force particles to stay increasingly close to the bottom of the wells. However, at any 
given e > 0, there is a positive probability that a particle switches from one well to the other 
in any given period of time. The rate at which this happens is governed by the local structure 
of H near ^ = ±1 and near ^ = 0, and becomes very small — of order e"^ exp(— 1/e), as we 
shall see below. 

In the limit e = 0, the behaviour of particles in the ^-direction is no longer recognizable 
as diffusional in nature. In the ^-direction a particle can only be in one of two states ^ = ±1, 
which we therefore interpret as the A and B states. Of the diffusional movement in the S,- 
direction only a jump process remains, in which a particle at ^ = — 1 jumps with a certain 
rate to position ^ = 1, or vice versa. 

1.5 Spatiochemical rescaling 

Since the jumping (chemical reaction) rate at finite e > is of order exp(— 1/e), the 
limiting reaction rate will be zero unless we rescale the system appropriately. This requires 
us to speed up time by a factor of £exp(l/e). At the same time, the diffusion rate in the 
x-direction remains of order 1 as e ^ 0, and the rescaling should preserve this. In order to 
obtain a limit in which both diffusion in x and chemical reaction in ^ enter at rates that are 
of order 1, we use the freedom of choosing the parameter r that we introduced above. 
We therefore choose r equal to 

Te := eexp(l/£), (5) 

and we then find the differential equation 

dtpe- ^^Pe-T,di_{d^p, + \ped(_H)=0 in ^'(73 X (0,00)), (6) 

which clearly highlights the different treatment of x and ^: the diffusion in x is independent 
of Te while the diffusion and convection in the ^-variable are accelerated by a factor r^. 

1.6 Switching to the density variable 

As is already suggested by the behaviour of the invariant measure 7^, the solution p^ will 
become strongly concentrated at the extremities {±1} of the ^-domain (—1,1). This is the 
reason why it is useful to interpret p^ as a family Pe(t, •) of time- dependent measures, instead 
of functions. It turns out that the densities u^{t, •) 

MV) :='^ 

d7£ 

of /5e(f, ■) with respect to 7^ also play a crucial role and it is often convenient to have both 
representations at our disposal, freely switching between them. In terms of the variable 
equation ([6| becomes 

dtu, ~ A^u, - r,(a|5 - ^d^Hd^u,) = in (0, +00) x D, (7) 

supplemented with the boundary conditions 

d^Ue{t,x,±l) =0 for aU a; G n, \/^u^{t, x, ^) ■ n = on dfl x [-1,1], t >0. (8) 
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We choose an initial condition 

u,{0,x,0 =u°{x,0, for all (a;, G D, with p'l ^ u°j, € .^(D). (9) 

Let us briefly say something about the functional-analytic setting. It is well known (see 
e.g. [7]) that the operator := — — T^d^^ + {t^ / £)H' with Neumann boundary condi- 
tions Q has a self-adjoint realization in the space Hn :— L^(Z);7e). Therefore the weak form 
of equation ^ can be written as 

b^{dtu{t),v)+a^{u{t),v) = for ah « G "14, (10) 

where the bilinear forms and are defined by 

be : Hs X Hi: — » R, be{u,v) := / uvA'ye, 

J D 

and 

V, := W^-\D--f,) := [u G L'(D,7.) n W^^^{D) : J \V,,^u\^ dj, < +00} , 

Since 14 is densely and continuously imbedded in , standard results on variational evolution 
equations in an Hilbert triplet (see e.g. |13ll6]) and their regularizing effects show that a unique 
solution exists in C([0, 00); H^) f] C°°((0, 00); V4) for every initial datum u2 G H^. 

1.7 Main result I: weak convergence of and 

The following theorem is the first main result of this paper. It states that for every time 
t > the measures Pe{t) solutions of ^ weakly-* converge to a limiting measure p(t) in 
S^{D), whose density u{t) = is the solution of the limit system |TJ. Note that for a 

function u G L^{D,'y) the traces it* = «(-,±l) G L^{Q) are well defined (in fact, the map 
u I— » (m~,u*) is an isomorphism between L^{D,'y) and L^{Q, |Af2;R^)). 

We state our result in a general form, which holds even for signed measures in ^{D). 

Theorem 1. Let = u^-y^ G C*^ {[0, +00); .-^(D)) be the solution of ([^[[^ with initial 
datum p". If 

sup / |ji°|'d7, < +00 (11) 

e>0 Jd 

and 

pi weakly-* converges to = ii''7 — -^n ® "^-i + ^""^ ® '^+1 is e .[ 0, (12) 

then u''' G L^{D,'y), u'^'* G L^(r2), and for every t>0 the solution pe{t) weakly-* converge to 
p{t) = u{t)-f^ ^ti"(f)AQ®(5_i + ^M+(f)An®(5+i, (13) 
whose densities belong to C"([0, -l-oo); {yL))r\C^ {{Q, -\-oo); W^''^{Q.)) and solve the system 

dtu^ - A^M+ = k(u~ ~ it+) in SI X (0, +00) (14a) 
dtu~ — Aa;U^ — k{u^ ^ u^) in Q X (0, +co) (14b) 
u'^{0)=u°''^ mn. (14c) 



The positive constant k m ( 14a b) can he characterized as the asymptotic minimal transition 
cost 

k = -^\H"mH"{l) = limminjr, ((^'(0)' d7, : <y= G W^'''(-l, 1), V'(±l) = ±i}. 

(15) 
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Remark (The variational structure of the limit problem) . The "e — 0" limit problem ( 14a 



14cl admits the same variational formulation of the "e > 0" problem we introduced in 
Section 1.6 Recall that 7 is the measure defined in Q as the weak limit of 7^; we set 
H :— L%D, 7), and for every p = wy with u € H we set u^{x) :— u{x,±l) € L'^(f2, An). We 
define 

b{u,v) ~ J u{x,^)v{x,^)d-y{x,(,) ^ ^ J (u^ + u~ dXn- (16) 

Similarly, we set V := {u G H : G W^'^{i})}, which is continuously and densely imbedded 
in H, and 

a{u,v)~^f (^V^u+V.iJ+ + V^m"V:,w" +fc(?i+ -m")(«+ -11")) dAn. (17) 



Then the system ( 14 i,b,c) can be formulated as 

b{dtu{t),v) + a(u{t),v) =0 for every t > and u G (18) 



which has the same structure as ( 10 1. 

1.8 Main result II: a stronger convergence of Us 

Weak-* convergence in the sense of measures is a natural choice in order to describe the 
hmit of pe, since the densities and the limit density u = (it+,it^) are defined on different 
domains with respect to different reference measures. Nonetheless it is possible to consider a 
stronger convergence which better characterizes the limit, and to prove that it is satisfied by 
the solutions of our problem. 

This stronger notion is modeled on Hilbert spaces (or, more generally, on Banach spaces 
with a locally uniformly convex norm), where strong convergence is equivalent to weak con- 
vergence together with the convergence of the norms: 

Xn —> X Xn ^ X and llsnll ^ ||a;||. (19) 

In this spirit, the next result states that under the additional request of "strong" convergence 
of the initial data u^, we have "strong" convergence of the densities Ue; we refer to |17l I10| 
(see also [21 Sec. 5.4]) for further references in a measure-theoretic setting. 

Theorem 2. Let p^, pi be as in Theorem^ If moreover 



lim6,(M°,«^) = ^(w^M^), (20) 



then for every t > we have 



and 



limb^{u^{t),Ui:{t)) = b{u{t),u{t)) (21) 

eiO 

lima^{u^{t),u^{t)) = a{u{t) , u{t)) . (22) 

eiO 

Applying, e.g., [2l Theorem 5.4.4] we can immediately deduce the following result, which 
clarifies the strengthened form of convergence that we are considering here. This convergence 
is strong enough to allow us to pass to the limit in nonlinear functions of : 

Corollary 3. Under the same assumptions as in Theorem^we have 

lim / f{x,^,u,ix,^,t))d-f,{x,0^ f f{x,^,u{x,^,t))d-f{x,^) (23) 
^io Jd Jd 

= (^f{x,-l,u~{x,t)) + f{x,l,u'^{x,t))^dXsi{x) for every t > 0, 

where f : D x M. ^ M. is an arbitrary continuous function satisfying the quadratic growth 
condition 

\fix,^,r)\<A + Br'^ for every (x,^) eD, r G E 
for suitable nonnegative constants A, _B G R. 
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1.9 Structure of the proof 

Let us briefly explain the structure of the proof of Theorems [l] and [2] This will also clarify 
the term F-convergence in the title, and highlight the potential of the method for wider 
application. 

The analogy between ( |10[ ) and ( |18[ ) suggests to pass to the limit in these weak formulations, 
or even better, in their equivalent integrated forms 



b!;{us{t),Vs) + / as{us{t),Ve) dt = b{u^,V£), b{u{t),v) + / a{u(t) , v) dt = b{u\ v) . 
Jo Jo 



(24) 



Applying standard regularization estimates for the solutions to ( 10 1 and a weak coercivity 
property of be, it is not difficult to prove that Ue{t) "weakly" converges to u{t) for every t > 0, 
i.e. 

Pe(t) — Ue{t)'ye pif) = ^(i)7 weakly-* in ./£[D). 

The concept of weak convergence of densities that we are using here is thus the same as in 
Theorem 111 i.e. weak-* convergence of the corresponding measures in 



In order to pass to the limit in ( 24 1 the central property is the following weak-strong 
convergence principle: 

For every v £ V there exists G 14 with ^ u as e — > such that for every 
Ue ^ u 

be{ue,Ve) ^ b{u,v) and ae{ue,Ve) a{u,v). 

Note that the previous property implies in particular that recovery family converges 
"strongly" to v, according to the notion considered by Theorenij2l i.e. Ve v IS ^ v with 
both be{v!,,Ve) b{v,v) and ae(we,?;E) — > a{v,v). Corollary |6| shows that this weak-strong 
convergence property can be derived from F-convergence in the "weak" topology of the family 
of quadratic forms 

q^iu) := be(u,u) + Kae{u,u) to q"^ (u) :— b{u,u) + k a{u,u) for k > 0. (25) 

In order to formulate this property in the standard framework of F-convergence we will extend 
tte and &£ to lower semi-continuous quadratic functionals (possibly assuming the value -|-oo) 
in the space ./^{D), following the approach of jS) Chap. 11-13]. While the F-convergence of 
be is a direct consequence of the weak convergence of 7e to 7, the convergence of is more 
subtle. The convergence of and the structure of the limit depends critically on the choice 
of Te (defined in (|5|): as we show in Section 3.2 the scaling of in terms of e is chosen 



exactly such that the strength of the 'connection' between ^ = —1 and ^ = 1 is of order 0(1) 
as e — » 0. 

The link between F-convergence and stability of evolution problems of parabolic type is 
well known when 6^ = 6 is a fixed and coercive bilinear form (see, e.g., [H Chap. 3.9.2]) and 
can therefore be considered as the scalar product of the Hilbert space He = H. In this case 
the embedding of the problems in a bigger topological vector space (the role played by ./£[D) 
in our situation) is no more needed, and one can deal with the weak and strong topology 
of H, obtaining the following equivalent characterizations (see e.g. [5] Th. 3.16] and 8, Th. 
13.6]): 

1. Pointwise (strong) convergence in H of the solutions of the evolution problems; 

2. Pointwise convergence in H of the resolvents of the linear operators associated to the 
bilinear forms a^; 

3. Mosco-convergence in H of the quadratic forms associated to a^; 

4. F-convergence in the weak topology of H of the quadratic forms b + Ka^ to b + Ka for 
every «: > 0. 

In the present case, where be does depend on e, F-convergence of the extended quadratic 
forms be + Kae with respect to the weak-* topology of ./£[D) is thus a natural extension of the 
latter condition; Theorem[4]can be interpreted as essentially proving a slightly stronger version 
of this property. Starting from this F-convergence result, we will derive the convergence of 
the evolution problems by a simple and general argument, which we will present in Section |4] 
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1.10 Discussion 

The result of Theorem[l]is amongst other things a rigorous version of the result of Kramers [T^ 
that was mentioned in the introduction. It shows that the simple reaction-diffusion sys- 



tem ( 14 1 can indeed be viewed as an upscaled version of a diffusion problem in an augmented 
phase space; or, equivalently, as an upscaled version of the movement of a Brownian particle 
in the same augmented phase space. 

At the same time it generalizes the work of Kramers by adding the spatial dimension, 
resulting in a limit system which — for this choice of t^, see below for more on this choice — 
captures both reaction and diffusion effects. 

Measures versus densities. It is interesting to note the roles of the measures Pe,p and 
their densities Us,u with respect to 7^,7. The variational formulation of the equations are 
done in terms of the densities u^,u but the limit procedure is better understood in terms of the 
measures Pe, p, since a weak-* convergence is involved. This also allows for a unification of two 
problems with a different structure (a Fokker-Planck equation for and a reaction-diffusion 
system for the couple u~, u'^ .) 

Gradient flows. The weak formulation ( |10[ ) shows also that a solution Ue can be interpreted 
as a gradient flow of the quadratic energy |ae(ii, u) with respect to the L^{D,yi;) distance. 
Another gradient flow structure for the solutions of the same problem could be obtained by 
a different choice of energy functional and distance: for example, as proved in [11], Fokker- 
Planck equations like ([6| can be interpreted also as the gradient flow of the relative entropy 
functional 

in the space £P{D) of probability measures endowed with the so-called i'^-Wasserstein distance 
(see e.g. [2]). Other recent work 1 suggests that the Wasserstein setting can be the most 
natural for understanding diffusion as a limit of the motion of Brownian particles, but in this 
case it is not obvious how to interpret the limit system in the framework of gradient flows on 
probability measures, and how to obtain it in the limit as e — > 0. 

In a forthcoming paper we investigate a new distance for the limit problem, modeled on 
the reaction-diffusion term, and we study how the limit couple of energy and dissipation can 
be obtained as a F-limit. 

The choice of t^. In this paper the time scale is chosen to be equal to eexp(l/e), and it 
is a natural question to ask about the limit behaviour for different choices of . If the scaling 
is chosen differently — i.e. if t^£~^ exp(— 1/e) converges to or 00 — then completely different 
limit systems are obtained: 

• If Tie ^ eexp(l/e), then the reaction is not accelerated sufficiently as e — > 0, and the 



limit system will contain only diffusion (i.e. = in (14l); 

• If Tj: ^ eexp(l/e), on the other hand, then the reaction is made faster and faster as 
e — > 0, resulting in a limit system in which the chemical reaction A ^ B is in continuous 
equilibrium. Because of this, both A and B have the same concentration u, and u solves 
the diffusion problem 

dtu — Am, for a; G SI, t > 

■u.{0,x) = ^{u°-+{x) +u°''{x)) for xen. 

Note the instantaneous equilibration of the initial data in this system. 

While the scaling in terms of e of can not be chosen differently without obtaining 
structurally different limit systems, there is still a choice in the prefactor. For := fee^^^ 
with f > fixed, the prefactor f will appear in the definition ( |15[ ) of k. 

There is a also a modelling aspect to the choice of r. In this paper we use no knowledge 
about the value of r in the diffusion system at finite e; the choice r = Tg is motivated by 
the wish to have a limit system that contains both diffusive and reactive terms. If one has 
additional information about the mobility of the system in the x- and ^-directions, then the 
value of T will follow from this. 
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Equal rate constants. The assumption of equal depth of the two minima of H corresponds 
to the assumption (or, depending on one's point of view, the result) that the rate constant 



in ( 14 1 is the same for the two reactions A B and B A. The general case requires a 
slightly different choice for H, as follows. 

Let the original macroscopic equations for the evolution of A and B (in terms of densities 
that we also denote A and B) be 

dtA- l\A = k-B-k+A (27a) 
dtB - AB = k+A- k'B. (27b) 

Choose a fixed function Hq G C°°([-l,l]) such that Ho{±l) = and Ho{l) - Ho{~l) = 
logfc" — logfc^. We then construct the enthalpy by setting 

Hs := Ho H — H, 

£ 

where H is the same enthalpy function as above. The same proof as for the equal-well case 



then gives convergence of the finite-e problems to (271 



Equal diffusion constants. It is possible to change the setup such that the limiting system 
has different diffusion rate in A and B. We first write equation Q as 

dtp~divD,F, =0, 

where the mobility matrix G and the flux are given by 

By replacing the identity matrix block / in by a block of the form a(^) I the a;-directional 
diffusion can be modified as a function of ^. This translates into two different diffusion 
coefficients for A and B. 

The function H. The limit result of Theorem [l] shows that only a small amount of 
information about the function H propagates into the limit problem: specifically, the local 
second-order structure of H around the wells and around the mountain-pass point. 

One other aspect of the structure of H is hidden: the fact that we rescaled the ^ variable by 
a factor of ^/t^ can also be interpreted as a property of H, since the effective distance between 
the two wells, as measured against the intrinsic distance associated with the Brownian motion, 
is equal to 2^/7> after rescaling. 

We also assumed in this paper that H has only 'half wells, in the sense that H is defined 
on [—1, 1] instead of R. This was for practical convenience, and one can do essentially the 
same analysis for a function H that is defined on R. In this case one will regain a slightly 
different value of k, namely k = ^\H"{0)\H"{1)/2tv. (For this reason this is also the value 
found by Kramers |12l equation (17)]). 

Single particles versus multiple particles, and concentrations versus probabilities. The 
description of this paper of the system in terms of a probability measure p on D is the 
description of the probability of a single particle. This implies that the limit object {u~ ,u'^) 
should be interpreted as the density (with respect to 7) of a limiting probability measure, 
again describing a single particle. 

This is at odds with common continuum modelling philosophy, where the main objects are 
concentrations (mass or volume) that represent a large number of particles; in this philosophy 



the solution {u~ , u ) of (141 should be viewed as such a concentration, which is to say as the 
projection onto s-space of a joint probability distribution of a large number of particles. 

For the simple reaction A ^ B these two interpretations are actually equivalent. This 
arises from the fact that A —> B reaction events in each of the particles are independent of 
each other; therefore the joint distribution of a large number of particles factorizes into 
a product of A'^ copies of the distribution of a single particle. For the case of this paper, 
therefore, the distinction between these two views is not important. 

More general reactions. The remark above implies that the situation will be different for 
systems where reaction events cause differences in distributions between the particles, such 
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as the reaction A + B ^ C. This can be recognized as follows; a particle A that has just 
separated from a B particle (in a reaction event of the form C ^ A + B) has a position that 
is highly correlated with the corresponding B particle, while this is not the case for all the 
other A particles. Therefore the A particles will not have the same distribution. The best one 
can hope for is that in the limit of a large number of particles the distribution becomes the 
same in some weak way. This is one of the major obstacles in developing a similar connection 
as in this paper for more complex reaction equations. 

1.11 Plan of the paper 

One of the main difficulties in the proof of Theorem[l] namely the singular behaviour given by 
the concentration of the invariant measure 7^ onto the two lines at ^ = ±1, can be overcome 
by working in the underlying space of (signed or probability) measures in D. This point of 
view is introduced in Section [2| Section [s] contains the basic F-convergence results (Theorem 
[4| and the proof of Theorem [ijand of Theorem [5] The argument showing the link between F- 
convergence of the quadratic forms a^, and the convergence of the solutions to the evolution 
problems (see the comments in section 1.9 1 is presented in Section [4] in a general form, which 



can can be easily applied to other situations. 

2 Formulation of the evolution equations in measure 
spaces 

The Kramers-Smoluchowski equation 

We first summarize the functional framework introduced above. Let us denote by (■,•)£ the 
scalar product in M'' x M defined by 

{x,y)s ■- X ■ y + T^£_r), for every a; = (x, C), y = (y, ??) G K'* x M, (28) 

with the corresponding norm || ■ We introduced two Hilbert spaces 

H,~L\D,'y,) and V, = W^'\D,^,), 

and the bilinear forms 



be{u,v) :— / uvd^e for every u,v £ H^, (29) 

J D 

ai;{u,v) :— / (Vj^^u, Ve,^?;)^ d7£: for every u, w G \4, (30) 

J D 

with which ^ has the variational formulation 

bsidtUs,^) + a.s{us,v) = for every u G 14, t > 0; U!,{0,-) = u°. (31) 



The main technical difficulty in studying the limit behaviour of (31 1 as £ | consists of 



the e-dependence of the functional spaces H^, V^. Since for our approach it is crucial to work 



in a fixed ambient space, we embed the solutions of (31 1 in the space of finite Borel measures 



^{D) by associating to the measure :— u^'y^- We thus introduce the quadratic forms 

be{p) ■-be{u,u) ii p<^'y, andu= ^ e H., (32) 

d7£ 

Oe(p) := ae{u,u) if p < 7e and u = G Ve, (33) 

d7e 

trivially extended to +00 when p is not absolutely continuous with respect to 7^ or its density 
u does not belong to He or Ve respectively. Denoting by Dom(aE) and Dom(be) their proper 
domains, we still denote by ae{-,-) smd be{-,-) the corresponding bilinear forms defined on 



Dom(ae) and Dom(b£:) respectively. Setting pe := Uefe, cr := v^fe, (311 is equivalent to the 
integrated form 



be{pe{t),^)+ / ae{pe{r),(T)dr ^ be(p'l,a) for every (T G Dom(ae 
Jo 



(34) 
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We also recall the standard estimates 



for every t > 0, (35) 



ta,{p4t)) + 2 rb,{dtpe{r))dr= 0^(^,(0) dr for every t > 0, (36) 



^b4pe{t))+ta,{p,{t))+th,(dtpe{t)) < ^b,{p'^) for every t > 0. (37) 

Although versions of these expressions appear in various places, we were unable to find a 
reference that completely suits our purposes. We therefore briefly describe their proof, and 
we use the more conventional formulation in terms of the bilinear forms and and spaces 

and V^; note that be is an inner product for H^, and be + tte is an inner product for Ve. 

When uq is sufficiently smooth, standard results (e.g. |^ Chapter VII]) provide the exis- 
tence of a solution Ue G C([0,c>o); Ve) H C°°((0, oo); Ve), such that the functions t i—> ae{ue{t)) 
and t be{dtUe{t)) axe non-increasing; in addition, the solution operator (semigroup) St is a 
contraction in He- For this case all three expressions can be proved by differentiation. 

In order to extend them to all Ug G He, we note that for fixed t > the two norms on He 
given by (the square roots of) 



IbeiUe) 



and 



he{StU°e)+ f ae{SrU°e) 
^ Jo 



dr 



(38) 



are identical by (35 1 on a He-dense subset. If we approximate a general it^ £ He by smooth 
Ue^n, then the sequence Ue^„ is a Cauchy sequence with respect to both norms; by copying 
the proof of completeness of the space L^{0, oo; Ve) (see e.g. [SJ Th. IV. 8]) it follows that the 
integral in (38 1 converges. This allows us to pass to the limit in (35 1. The argument is similar 
for (371, when one writes the sum of (35 1 and (36 1 as 



lbe{ueit)) + tae{ne{t)) + 2 f rfe, (ftu, (r)) dr = J 6, (m^) . 



(39) 



Finally, (|37[) follows by (|39[) since r i— > be{dtUe{r)) is non-increasing. 



The reaction-diffusion limit 

We now adopt the same point of view to formulate the limit reaction-diffusion system in the 
setting of measures. Recall that for u E H := L^(D,7) we set u^{x) := u(z,±l), and thus 
we defined the function space 

V ■- {u (£ H : u"^ £ W''^(n)}, 

and the bilinear forms 



b{u,v) = ^f (u'^v'^+u V ^ dAn, 



a{u, v) := 



2 Jn ^ 



dXn. 



As before we now extend these definitions to arbitrary measures by 



b(p) := b{u,u) 
a{p) := a{u, u) 



■ r I dp 

it p J and u — — € H , 
d7 

if p <C 7 and u — dV , 
d7 



(40) 
(41) 

(42) 
(43) 



with corresponding bilinear forms b(-, ■) and o(-, ■); problem ( 14a b,c) can be reformulated as 

b{dtp(t), a) + a{p{t),a) — for every t > and a G Dom(a), 
or in the integral form 



b(p{t),a)+ / Q(p(r), cr) dr = b(p°, a) for every ct £ Dom(a) 
Jo 



(44) 



Since both problems (34 1 and (44 1 are embedded in the same measure space ^(D), we 
can study the convergence of the solution p^ of (341 as e | 0. 
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3 F-convergence result for the quadratic forms ag, be 



The aim of this section is to prove the following F-convergence result: 
Theorem 4. If p as e [ Q in ^{D) then 



liminf Ogfoe) > a(p), limini brips) > b(p). 



(45) 



For every p £ ^{D) such that a{p) + b{p) < +oo there exists a family p^ £ ./£[D) weakly-* 
converging to p such that 



\imae{pe) = a{p), lim be(pe) = b(p). 



(46) 



Note that ^{D) endowed with the weak-* topology is the dual of a separable Banach 
space, and therefore the sequential definition of F-convergence coincides with the topological 
definition [8] Proposition 8.1 and Theorem 8.10]; consequently Theorem |4] implies the F- 
convergence of the families and be. Theorem |4] actually states a stronger result, since the 
recovery sequence can be chosen to be the same for and b^. This joint F-convergence of 
the families and be is nearly equivalent with F-convergence of combined quadratic forms: 

Lemma 5. Theorem^ implies the 

V{^{D))- convergence of cC{p) be{p) + na^ip) to c{'^{p) b{p) + Ka{p) (47) 
for each k > 0. 



Conversely, if we assume (47 1, then (461 holds, and (451 follows under the additional 
assumption 

limsupae(pe) + be(pe) = C < +00. (48) 



Proof. The first part of the Lemma is immediate. For the second part, suppose that p^ 
and satisfies (481; the F-liminf inequality for q? yields 

liminf be(pe) > liminf q^ipe) — Ck > q''(p) — Ck = b(/5) + k (a(p) — C) for every k > 0, 



and therefore the second inequality of (451 follows by letting k j 0. A similar argument yields 
the first inequality of ( 45 I . 

Concerning (461, F-convergence of to q^ yields a recovery family p^ ^ p such that 

Xvcao.e'yPe) + be(pe) = a(p) + b(p) < +00; 



In particular o.e(^pe^^be{pe) is uniformly bounded, so that (45 1 yields the separate convergence 
(il. □ 



One of the most useful consequences of (47 1 is contained in the next result (see e.g. |16l 
Lemma 3.6]). 



Corollary 6 (Weak-strong convergence). Assume that (47 1 holds for every k > and let 
Pe,(Je G ./#(!)) be two families weakly converging to p,a as e [ and satisfying the uniform 
bound (48 1, i.e. 



limsupae(pe) + be[ps) < +00, limsupae(aE) + b^ioe) < +00, 

ej.0 £iO 

so that p, a belong to the domains of the bilinear form a and b . We have 



(49) 



livaaeioA = a(cr) 
lim beioe) ~ b(a) 



\uaae{pe,crs) = o(p, a) 

lim be (pejCTe) = b(p, Cr). 
ei-O 



(50) 
(51) 
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Proof. We reproduce here the proof of ^16] in the case of the quadratic forms (501. Note 
that by (49 1 and Lemma [s] we can assume that and satisfy (451. For every positive 
scalar r > we have 



Taking the inferior hmit as e | and recalhng (451 we get for A := hmsup^^Q ae{pe) 
Um^inf 2ae(pe, os) > a{rp + r^^a) — r^A — r~^a(a) = 2a(p, a) + (o(p) — A). 



Since r > is arbitrary and A is finite by ( |49[ ) we obtain hminf^jo Oe(pe, erg) > a{p,a) and 
inverting the sign of a we get (|50[). □ 



We spht the proof of Theorem [4] in various steps. 

3.1 Estimates near x { — 1, 1}. 

Lemma 7. If Pe = Us^e satisfies the uniform bound aei^ps) < C < +cx3 for every e > 0, then 
for every 5 G (0, 1) 

d^u^ -^0 in L^(Q x tug), as e — > 0, (52) 
where ujs := (—1, —S) U {S, 1) . 

Proof. We observe that 

Te / (95U£)^d7£: < Oeips) <C<00. 
J D 

If fti = sup H{^) < 1, then inf e'"'-^'>^^ = e~''■"'^ and we find 

^2J At ^ r^'^^ ^ r^^-^ 
(d^Mj:) d3;d^<G — =G — e = . 



Taking the hmit as e ^ we obtain (52 1. □ 

Lemma 8 (Convergence of traces). Let us suppose that p^ = Us:"/^ p = wy with ae{ps) < 
C < +00 and let uf{x) be the traces of at £^ — ±1. Then as e I 

u^ strongly in L^{Q), (53) 

where are the functions given by (|13[) . 



Proof. Let us consider, e.g., the case of . Let us fix 5 £ (0, 1); by (521 and standard trace 
resuhs in W^'^{~1, —1 + 5) we know that 



r r-l+S 

lira / ijj^{x)dL =0 where u>^{x) :— sup \u£{x, £,) — u~ {x)\ <S \d^Us(x,£^)\ d^. 
<^ioJn -i<5<-i+i J_i 

(54) 

Let us fix a function tp e C°{n) and a function tp G C°[-l, 1] with < ^ < 1, tpi-l) = 1, 
snpptp C [—1, —1 + S]; we set 



' ^iOdUO, M^) f u,{x,0^{Od%{0, 

1 ./-I 



where 7^ is the measure defined in ([3|. Note that 



lim J, = (^,7) = i^(-l) + ^^{1) = 



Since Pe weakly converge to p we know that 

lim / (f){x)u^{x)dXn ^limJ^^ / (l>{x)^p{£,)us{x, £,) dj^ix,^) = / <;!>(x)u" (x) dAf 
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so that converges to u in the duahty with bounded continuous functions. On the other 
hand, 

/ \V^u,{x)fdXn < J-^ f r jV,«,(a;,C)|V(C)d7(C)dAo(x) < Jr'o,(p,) < 2C 
Jn JnJ-i 

so that -Se — > u~ in L'^{Q) by Rellich compactness theorem. 
On the other hand, thanks to (541, we have 

hm / Ljr(3;) - ile(a;)j dXn{x) = lira J^'^ f \ f 'tp{^){ui:{x, ^) - u~ (x)) d^s{^)\ dAn(a;) 

<hm / ^(CK2(x)d7.(x,e) =0, 



which yields ( 53 1 . □ 

Remark. A completely analogous argument shows that if p,. satisfies a W^'^{D; 7E)-uniform 
bound 

/ ||V.,eii,||,d7e(x,0 < C < +00 (55) 

J D 

instead of a^ips) < C, then — > ii^ in L^{Q). 

3.2 Asymptotics for the minimal transition cost. 

Given {ip^,(p'^) £ let us set 

K,{^- ,^+) ~ mm[T, J'\^'iO)^ d^, : ^ € W''\-l,l), ^(±1) = <^±} (56) 

It is immediate to check that is a quadratic form depending only on — <^^, i.e. 

K,{^-,ip+) = k,{^+ ->fi'f, k,=K,{-l/2,l/2). (57) 

We call 7^{(fi~ , (f^) the solution of the minimum problem ( |56[ ): it admits the simple repre- 
sentation 

T,(v5-,</p+) = ^(^"+^+) + (<^+-¥'-)<^e (58) 
where (p^ = 7^{~l/2, 1/2). We also set 

Q.(^-,^"') (T.(^-,^+))'d7. = ^((^-)' + (¥>+)') + (g. - !)(¥'+ - ^-)' (59) 

where ^ 

:= I ^ 10.(01' d7,(C) = Q.(-l/2, 1/2). (60) 

Lemma 9. We have 

k ^-H"{0)H"{1) 

'm''^ = -2 = — — ' (6^) 



and 



limge = ^ so that limQeii^} ,V'^) = ^(</^ f + (62) 



Proof. solves the Euler equation 

(e-"«)/^,^Ue))' = on (-1,1), 0,(±1) = ±|. (63) 



We can compute an explicit solution of ( 63 1 by integration: 



Jo 



drj. 
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Define :— J^^ e^(€)/= d^. The boundary conditions for ^ = ±1 give 

C' = 0, C J e"^^'^'' = C/, = 1 

It follows that 

Jo 

and 

We compute, using Laplace's method: 



^1^6^/^(1 + 0(1)) and = J-|^(l + o(l)), ase^O, 



thus obtaining ([61|. Since 

^: = I-'e"^' ^5o in ^'(-1,1) 

and H is even, we have 

<^.(C) = 4"' ^^""^"^/^ d,) -> ^sign(e) 

uniformly on each compact subset of [—1, 1] not containing 0. Since the range of (fig belongs 
to [-1/2, 1/2] and % + we obtain (|62|. □ 

3.3 End of the proof of Theorem |4| 

The second limit of ( |45[ ) follows by general lower semicontinuity results on integral functionals 
of measures, see e.g. [2) Lemma 9.4.3]. 

Concerning the first "lim inf " inequality, we split the the quadratic form in the sum of 
two parts. 



We choose a smooth cutoff function rj : [— 1, 1] ^ [0, 1] such that -q (— 1) = 1 and supp(»7 ) C 
[— 1, —1/2] and the symmetric one r)^{^) := vi^O- We also set 

u:{x):^ y\-(eK(a;,C)d7.(0, utix):^ j' {Ou,{x,0 d%{0, (65) 
and it is easy to check that 

uf ^ ^u"^ in (66) 

We also set 6^ := f^^V^ ( = jl^V' (0 <i%iO) ■. observing that 6^ 1/2. We then 

have by Jensen inequality 

al{Pe)> I /\?7-(C) + ^?+(C))|V.«.|'d7,(C)dAn > / \V \' + \V ut \^ dXn 
and, passing to the limit, 

liminf Oe(p£:) > ^ / I Vit^ P + I Vu^ P dAn 
eio ' - 2 in ' 



Let us now consider the behaviour of o^: applying ([56f and | |57[ ) we get 

ai{pe) = J (je J {d^Us(x,^)f d^siO^ d\n> J ks{u^ {x) - uf (x)f dXn 
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so that by (61 1 and (531 we obtain 



k z 
liminf a^(/9e) > - / (u~ (x) — u'^ (x))^ d\n- (67) 



In order to prove the "hm sup" inequahty {46\ we fix p — u-y with u in the domain of the 



quadratic forms a and b so that u = u{-, ±1) belong to W ' {Q), and we set = Ueje where 



Ue{x,-) = 7s{u {x),u~^{x)) as in (581. We easily have by (62i and the Lebesgue dominated 
convergence theorem 

limbe:(pE) = hm / Qe(u"(a;), u+(a;)) dAn = / (]r\u^ (x)]'^ + hu'^ {x)\^) dXn = b(p). 

Similarly, since for every j = 1, • • • ,d and almost every x £ Q 

dxjUe{x,S) = 7{dxjU~ {x),dx.u'^), 

we have 

„ d 

hmOe(p£) = hin / (y^^Qe{da:-u~{x),dx.u^{x)) + K^[u^{x),u'^{x))) dXn = 
^ ^ Jn j^i 

= 1^ {l\^u-ix)\' + l\\/u+{x)f + '^{u-{x) - u+{x))) dXn = n(p). 

4 From F-convergence to convergence of the evolu- 
tion problems: proof of Theorems 1 and 2. 

Having at our disposal the F-convergence resuh of Theorem |4] and its Corollary |6] it is not 
difficult to pass to the limit in the integrated equation (34 1. 

Let us first notice that the quadratic forms be satisfy a uniform coercivity condition: 

Lemma 10 (Uniform coercivity of b^). Every family of measures p^ £ ^{D), e > satisfying 

limsupb£(pe) < +00 (68) 

e>0 

is bounded in ^{D) and admits a weakly-* converging subsequence. 

Proof. It follows immediately by the fact that 7^ is a probability measure and therefore 



( 68 1 thus implies that the total mass of p^ is uniformly bounded and we can apply the relative 



weak-* compactness of bounded sets in dual Banach spaces. □ 

The proof of Theorems [l] and [2] is a consequence of the following general result: 

Theorem 11 (Convergence of evolution problems). Let us consider weakly-* lower-semicontt- 
nuous, nonnegattve and extended-valued quadratic forms a^, be, a, b defined on .^(D) and let 
us suppose that 

1) Non degeneracy of the limit forms: b is non degenerate (i.e. b{p) — p = 0) 
and Don (a) is dense in Dom(b) with respect to the norm-convergence induced by b. 

2) Uniform coercivity: b^ satisfy the coercivity property stated in the previous Lemma 

m 

3) Joint F-convergence: :— be + kOe satisfy the joint F-convergence property (|47[) 



r(.^(D))-limq^ = q"" = b Ka for every k>0. (69) 
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Let Pe{t), t > 0, be the solution of the evolution problem (34l starting from p^' G Dom(b£). 
// 

^ m ./£{D) as e [Q wtth limsupbe(pe) < +00 (70) 

then Pe(i) ^ p(i) in ^{D) as e [ Q for every t > G and pit) is the solution of the limit 
evolution problem ( 44 1 . 



// moreover Mvn^yo be(p?) — b(po) then 

\\rahi:{pe(t)) = \}[p(t)), limae(pe(t)) = o(p(t)) for every t > 0. 



(71) 



Proof. Let us first note that by (351 and the coercivity property of be the mass of Pe(t) is 
bounded uniformly in t. Moreover, (37 1 and the coercivity property show that dtps is a finite 
measure whose total mass is uniformly bounded in each bounded interval [fo,ii] C (0, +cx3). 
By the Arzela-Ascoli theorem we can extract a subsequence p^^ such that p^^ (t) pit) for 
every t>Q. The estimates ( 37 1 and ( 45 1 show that for every t > 0, p(t) belongs to the domain 
of the quadratic forms and b, and satisfies a similar estimate 

\b{p{t)) + ta{p{t))^t^b{dtp{t)) < i^liminf b(p^) < +00. (72) 

2 2 eiO 



Let a £ ^{D) be an arbitrary element of the domains of a and b; by (47 1 we can find 
a family (actually a family , but we suppress the subscript n) weakly converging to a 
such that (461 holds. By (341 we have 



be(pe(t),(T£) + / Oe(pe(r),(Je)dr = be(p°,0-£) 



(73) 

and ( |37[ ) with the Schwarz inequality yields the uniform bound 

jae(pe(t),ae)j < t-^'^b,{plf'\,{a,f^'' < Ct~^'^ 
where C is independent of e; we can therefore pass to the limit in (73 1 by Corollary |6] to find 

b{p{t),a)+ f ne(p(r),cr)dr = b(po,(j), 



so that p is a solution of the limit equation. Since the limit is uniquely identified by the 
non-degeneracy and density condition 1), we conclude that the whole family p^ converges to 
p as £ I 0. In particular p satisfies the identity 



Jb(p(i)) + / o(p(r)) dr = Jb(p") for every t > 0. 



(74) 



This concludes the proof of ( 70 1 (and of Theorem . 

In order to prove ( |71[ ) (and Theorem [2| we note that by ( |35| l and ( |74[ ) we easily get 

limsup Jb,(pe(t))+ / a,{pe{r))dr<h{p{t))+f a{p{r))dr. 



The lower-semicontinuity property ( 45 1 and Fatou's Lemma yield 

o(p(r)) dr for every t > 0. (75) 

□ 



Iimb,(pe(t)) = b(p(t)), lim/ aApe{r))dr 



Applying the same argument to (36 1 and its "e = 0" analogue we conclude that as{pe{t)) 
a{p{t)) for every t > 0. 



Remark (More general ambient spaces). The particular structure of ^{D) did not play any 
role in the previous argument, so that the validity of the above result can be easily extended 
to general topological vector spaces (e.g. dual of separable Banach spaces with their weak-* 
topology), once the coercivity condition of Lemma 10 is satisfied. 
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